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We study the generation of parallel electric fields by virtue of propagation of ion cyclotron waves 
(with frequency 0.3 cu ci ) in the plasma with a transverse density inhomogeneity. Using two-fluid, 
cold plasma linearised equations, we show for the first time that En generation can be understood by 
an analytic equation that couples E\\ to the transverse electric field of the driving ion cyclotron wave. 
We prove that the minimal model required to reproduce previous kinetic results on E\\ generation 
is the two-fluid, cold plasma approximation in the linear regime. In this simplified model, the 
generated E\\ amplitude e.g. for plausible solar coronal parameters attains values of w 90 Vm _1 for 
the mass ratio mi/m e = 262, within a time corresponding to 3 periods of the driving ion cyclotron 
wave. By considering the numerical solutions we also show that the cause of Eu generation is 
electron and ion flow separation (which is not the same as electrostatic charge separation) induced 
by the transverse density inhomogeneity. The model also correctly reproduces the previous kinetic 
results in that only electrons are accelerated (along the background magnetic field), while ions do 
not accelerate substantially. We also investigate how En generation is affected by the mass ratio 
and found that amplitude attained by Eu decreases linearly as inverse of the mass ratio mi/m e , i.e. 
En oc 1/rru. This result contradicts to the earlier suggestion by Genot et al (1999, 2004) that the 
cause of En generation is the polarisation drift of the driving wave, which scales as oc m,i. Also, for 
realistic mass ratio of mj/m e = 1836 our empirical scaling law is producing En = 14 Vm -1 (for solar 
coronal parameters). Increase in mass ratio does not have any effect on final parallel (magnetic field 
aligned) speed attained by electrons. However, parallel ion velocity decreases linearly with inverse 
of the mass ratio mi/m e , i.e. parallel velocity ratio of electrons and ions scales directly as mi/m e . 
These results can be interpreted as following: (i) ion dynamics plays no role in the Eu generation; 
(ii) decrease in the generated parallel electric field amplitude with the increase of the mass ratio 
nn/me. is caused by the fact that «<j = 0.3oj C i cx 1/mj is decreasing, and hence the electron fluid 
can effectively "short-circuit" (recombine with) the slowly oscillating ions, hence producing smaller 
E\\ which also scales exactly as l/roj. 

PACS numbers: 52.20.-j,52.25.Xz,52.30.Ex,52.35.-g,96.60.-j,96.60.Hv 



I. INTRODUCTION AND MOTIVATION 

The generation of parallel electric fields in inhomoge- 
neous plasmas is a generic topic, which is of interest in 
a variety of plasma phenomena such as particle accelera- 
tion in Solar and stellar flares [l[ , auroral acceleration re- 
gion and current sheets in the Earth magnetosphere (see 
refs. in Q), laboratory plasma reconnection experiments 
[HQ and many more. In situ and remote observations of 
accelerated particles often show parallel electric fields in 
localised double layers, charge holes or U-shaped voltage 
drops. 

In many astrophysical plasmas, an adequate form of 
description of large-scale, bulk dynamics is provided by 
Magnetohydrodynamics (MHD). However, MHD cannot 
provide proper description of some fundamental ques- 
tions such as dissipation (which necessarily occurs at 
small-scales) and particle acceleration, unless the con- 
cept of somewhat uncertain from the fundamental point 
of view anomalous resistivity is invoked. The particle ac- 
celeration is of a considerable importance e.g. for Solar 
flares where the accelerated particles gain 50-80% of the 
energy released during this process. On one hand, ob- 
servable dynamics e.g. (i) MHD waves in the case Solar 



plasmas; (ii) jets and accretion disks, in the case of stel- 
lar or compact objects or centres of Galaxies; and (iii) 
MHD waves in Tokamak spectroscopic studies; are well 
described by MHD theory. On the other hand, small- 
scale processes such as dissipation and particle accelera- 
tion are not observable directly. This creates controversy 
around issues such as the coronal heating problem (as to 
why the Solar corona is 200 times hotter than underly- 
ing photosphere); anomalous resistivity which manifests 
itself in an unusually fast damping of kink oscillations 
of solar coronal loops; or anomalous viscosity (problem 
of getting rid of angular momentum) in accretion disks. 
This dichotomy is schematically sketched in Fig. 1. Here 
energy cascade from the large scales to small scales is de- 
picted as either 1// = A; -1 the white noise spectrum (in 
the case of waves) or some form of turbulence spectrum 
(with some power law of k~ a dependent on a particular 
turbulence model). 



When MHD is used for the description of plasmas, the 
electric field is totally eliminated from the consideration. 
On one hand, this has a good justification due to the 
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FIG. 1: The sketch of typical power law spectrum of fluctua- 
tions along with their observability criteria. 
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i.e. for non-relativistic (V <C c) plasmas the ratio of 

the displacement and (V x B) currents is much smaller 
than unity. Note that in the Eq.(l) spatial and time 
derivatives were approximated by: d/dx rs 1/L and 
d/dt w 1/r, where L and T are typical spatial and 
temporal scales of the system; and the ideal MHD limit 
(E = —V x B/c) was used. Thus, by neglecting the 
displacement current the electric field is totally excluded 
from the consideration. On one hand, this assumption 
may well be valid for the large scales. On the other 
hand, when the small scales are considered, the electric 
field, which appears (as we will show below) because of 
the electron and ion flow separation (which is impossible 
to treat correctly in single fluid MHD) starts to play far 
more important role than previously thought. 

Authors of Ref. Q pointed out that previous studies 
of the £j| generation, based on the balance of the differ- 
ent terms in the generalised Ohm's law, were not properly 
addressing the issue. In essence their main argument was 
that in such approach the generalised Ohm's law merely 
states the Newton's second law F = ma, whilst obscuring 
the true source of the parallel electric field generation. It 
was suggested that the source of En is the parallel dis- 
placement current. As stated above, this term is usually 
ignored, however in the regions of low density, for a cer- 
tain (V x i?)|| , the plasma is too dilute to carry significant 
J|i and thus (l/c)dEn/dt becomes important [2J. One of 
the main conclusions that immediately follows is that the 
signatures of the generated En in space plasmas should 
be correlated with low plasma density. 

Yet another series of works exist, which investigate the 
generation of parallel electric fields by virtue of propa- 



gation of Alfven waves (or more precisely ion-cyclotron 
waves, see below) in the plasma with a transverse den- 
sity inhomogeneity [1, @, 0, 1, ©, EH EH- To this da y 
the true cause of the generation of En in these studies 
eluded determination. Authors of Ref. [5| considered the 
case of both transverse and longitudinal density inhomo- 
geneity, applicable to the stratified Earth magnetosphere. 
They demonstrated that E\\ is generated in the regions of 
transverse density gradients, and presented an analytical 
model in which the En and E± are coupled via longitu- 
dinal density gradient (see Eq.(6) from Ref. Q). Subse- 
quently, detailed numerical study of long term evolution 
of the system was presented, including the generation of 
E\\ @| . However, in Ref.0] only the case of transverse 
density inhomogeneity was considered, while theoretical 
explanation was still based on Ref. [5| . This seems incor- 
rect because the latter reference attributes E\\ and E± 
coupling to the longitudinal inhomogeneity, which is ab- 
sent in Ref. [6j . In brief, these two works suggest that the 
Alfven wave propagation on sharp density gradients leads 
to the formation of a significant parallel electric field. It 
results from an electric charge separation generated on 
the density gradients by the polarisation drift associated 
with the time varying Alfven wave electric field [6j . Their 
approach involved substituting ion polarisation drift cur- 
rent (electron one was omitted because of its proportion- 
ality to the particle mass) j± — (m^i/ B 2 )dE±/dt into 
the Maxwell equations, which with the aid of the conser- 
vation laws yielded the equation for En and E± coupling 
H . Unaware of these works authors of Refs. consid- 
ered similar physical system with the increased density in 
the middle of the domain (mimicking) solar coronal loop, 
as opposed to Earth magnetospheric density cavity case 
studied in Refs. 0,0]. Similar effect of En generation was 
found because of the existence of density gradients in the 
system. Later a comment paper was published 0] , which 
detailed similarities and differences of the two series of 
works. 

It should be noted in passing that at that time we 
came to the realisation that electron acceleration seen in 
both series of works 0, 0, 0, § is a non-resonant wave- 
particle interaction effect. In Refs. @, @] the electron 
thermal speed was v t h,e — 0.1c while the Alfven speed 
in the strongest density gradient regions was va — 0.16c; 
this unfortunate coincidence led us to the conclusion that 
the electron acceleration by parallel electric fields was af- 
fected by the Landau resonance with the phase-mixed 
Alfven wave. In Refs. 0, @ the electron thermal speed 
was Vth.e = 0.1c while the Alfven speed was va = 0.4c be- 
cause they considered a more strongly magnetised plasma 
applicable to Earth magnetospheric conditions. Based 
on this observation, Refs. [l(J [ll[ explored the possibil- 
ity of En generation in the MHD description in the solar 
coronal heating problem context. Although, in the lat- 
ter approach, the heating aspect seems certain (because 
the fast magnetosonic waves, which are generated by the 
interaction of weakly non-linear Alfven wave with the 
transverse density inhomogeneity, dissipate on the bulk 
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Braginkii resistivity), the issue whether such E\\ can ac- 
celerate particles is less clear [Til ]. 

II. THE MODEL AND RESULTS 

The above discussion demonstrates that the issue of 
true cause of Eu generation when an Alfven wave moves 
in the transversely inhomogeneous plasma eluded identi- 
fication. In this work we present a minimal model which 
can explain Eh generation in mathematically and phys- 
ically rigorous manner. We start from two-fluid, cold 
jignoring thermal pressure) plasma linearised equations 

d t Ve = -(e/m e ) (E + V e x £ /c) , (2) 
d t Vi = +(e/n»i) (E + V % x S /c) , (3) 

d t B = -cV x E, (4) 

d t E = cV xB-±7rne(Vi-V e ). (5) 

Hereafter subscripts under d denote partial deriva- 
tive with respect to that subscript. Uniform, 
background magnetic field, Bq is in z-direction. 
Density profile is specified as a ramp, n(x) — 
n (1 + 3exp [-[(x - 1005) / (208)] 6 ]) in which the cen- 
tral region (along x-direction, i.e. across z), is smoothly 
enhanced by a factor of 4, and there are the strongest 
density gradients having a width of about 206 around 
the points x — 81(5 and x = 119(5. Here 5 = c/oj pe is 
the (electron) skin depth, which is a unit of grid in our 
numerical simulation. We use 2.5D description meaning 
that we keep all three, x, y, z components of all vectors, 
however spatial derivatives d/dy = 0. The above nor- 
malised plasma number density and Alfven speed profiles 
are shown in Fig. (2). 

In order to derive the equation that describes Eu — 
E z generation, we write Eqs.(2)-(5) in x,y,z component 
form. Omitting details of the calculation we present the 
final result: 

(4 - c 2 8 2 xx + w£ + £,| = -c 2 d 2 zx E x . (6) 

Also, a similar calculation enables us to obtain the equa- 
tion describing the dynamics of driving transverse electric 
field E x of an ion cyclotron wave: 

{d 2 t -c 2 d 2 zz + uj 2 pi +J 2 pe )E x = 

~ <? d lx E \\ ~ uli(mi/e)u ci Viy - uj 2 pe {m e / e)uj ce V ey . (7) 

In the considered problem E x and B y are both compo- 
nents of Alfven (ion cyclotron) wave, so these can be used 
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FIG. 2: Dimensionless number density, open squares, and 
Alfven speed, solid line, profiles across the uniform unper- 
turbed magnetic field (i.e. along x-coordinate) which is used 
as an equilibrium configuration in our model of a footpoint of 
a solar coronal loop or a polar region plume. 

interchangeably. Note that Eq.(7) also describes the feed- 
back of the generated Eh on the driving transverse elec- 
tric field E x (see the first term on the right-hand-side). 
Here the notation is standard: u pe — \J Aim^e 2 / m e and 
ujpi = Aim^e 2 / rrii are electron and ion plasma frequen- 
cies; u> c (e,i) — e-Bo/( m (e,i)c) are respective cyclotron fre- 
quencies. 

It is interesting to note that Eqs.(6) and (7) can be 
also obtained from the dielectric permeability tensor of 
cold, magnetised plasma (e.g. chapter 4.9 in Ref.[l2|). 
For example, Eq.(6) can be directly obtained from the 
classical equation for the electric field perturbation, E\, 
in the case of Eq = and Bq = Bqz 

V x V x Ei = (aj 2 /c 2 )eE u (8) 

where e is the dielectric permeability tensor of cold, 
magnetised plasma. In effect, Eq.(6) can be obtained 
from the z-component of Eq.(8) and putting in e zz — 
1 - uj 2 Ju 2 - uJ 2 }i /uJ 2 . 

In order to solve Eqs.(2)-(5) numerically we use the 
following normalisation: t — tu>~,}, V x ^ VtZ = V x ^ VtZ c, 

Bx.y,z — E %,y ,z (jYleCLdpe/&) — E x ,y. z E() : B x .y jZ — 

B x ,y,zBo 7 and {x,y,z) — c/uj pe (x,y,z). In what follows 
we omit tilde on the dimensionless quantities. The (x, z) 
simulation 2D box size is 2005 x 2500(5. Since we fix back- 
ground plasma number density at 10 9 cm -3 (typical value 
for the solar corona), aj pe is then 1.784 x 10 9 rad s _1 and 
the simulation box size is 33.6 m in x- and 420.5 m in z- 
direction. Bo was fixed at 101.5 Gauss (typical value for 
the solar corona), which gives oj ce /uj pe = 1. rrii/m e ratio 
was varied as: 45.9, 91.8, 183.6, 262.286 (realistic one 
is 1836). These values correspond to 1/40, 1/20, 1/10 
and 1/7-th of the realistic value respectively. This yields 
respectively: u! C i/uj p i = Bq / {c^J Aiminii) — Va/c = 
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l/y/rru/me = 0.148, 0.104, 0.074 and 0.062 for x < 70 
and x > 130 (realistic to ci /uj p i = V/Jc is 0.023). Here pa- 
rameters are similar to e.g. Refs.UHj], except for far more 
realistic mass ratios. Note that the simulation parame- 
ters are still somewhat artificial. Full kinetic, Particle- 
In-Cell (PIC) simulations employed in Refs.0, Q or in 
gyro-kinetic approach which uses guiding centre approx- 
imation for electrons, whilst retaining ion particle- like 
dynamics @, [f| are computationally challenging. Thus, 
in those studies rather modest mass ratios e.g. 16 were 
used. Note also, that since here we do not need to re- 
solve electron thermal motions as we are only studying 
electromagnetic part of the problem (E\\ generation) our 
unit of spatial grid size is 5 = c/uj pe , the (electron) skin 
depth. While in full kinetic, PIC simulation @, @] the 
unit of grid has to be A = v t h,e/u pe . Since in a PIC sim- 
ulation typically v t h, e /c = 0.1, in the present, two-fluid 
approach an equivalent to PIC numerical simulation re- 
quires (S/A) 2 = {c/v t h,e) 2 — 10 2 less grid points, thus 
it can be 100 times faster. For comparison a single run 
for mass ration 16 in Refs.0, H| takes about 8 days on 
parallel, 32 dual-core 2.4 GHz Xeon processors, similar 
run with mass ratio of 262 would have taken 4 months. 
The numerical run presented here for the mass ratio of 
262 takes 4 days with only one processor. 

We solve relativistic version of Eqs.(2)-(5) numerically 
with a specially developed and tested FORTRAN 90 code 
which uses 4-th order centred spatial derivatives and 4- 
th order Runge-Kutta time marching. Although Alfven 
speeds considered are at most w 15 % of the speed of light 
for mi/me — 45.9, relativistic effects were included. The 
simplest option becomes available in the linear regime. In 
ref. [l2T |. appendix I, paragraph 5, it was shown that the 
relativistic equation of motion of a particle with charge 
q and the rest mass mo can be written as 



q 



TV 

dt 7777,0 



E 



V x B V(V ■ E) 



(9) 



where 7 = (1 — V 2 / c 2 )^ 1 / 2 . As can be seen from the 
latter equation, in the linear regime, it coincides with 
either Eq.(2) or (3) after substituting m e ^ — > 7 e ,i"ie,i> 
where 7^ = (1 — V^/c 2 ) -1 / 2 . Naturally, such simplified 
approach is only valid when there are no flows in the un- 
perturbed state Vo = 0. As can be seen below, largest at- 
tained velocities in the simulation are those of electrons, 
and these do not exceed 3 % of speed of light. Thus, rel- 
ativistic corrections play only a minor role. It should be 
noted, however we still retain the displacement current 
in Eq.(5). Note, also that the gradients in the code are 
resolved numerically to an appropriate precision (20 grid 
points (open squares) across each gradient in Fig. (2)). 

Initially all perturbations are set to zero, 
and we start driving the z = 1 cell with 
the transverse magnetic fields of the form 
B y = -0.05 sin(w d i) (1.0- exp[-(</(3.125w ( :/)) 2 ]) and 
B x = -0.05cos(w d i) (l.0-exp[-(t/(3.125w ( T i 1 )) 2 ]). 
As in Refs.0, Q, we fixed aid at 0.3 oj C i (to 



avoid ion-cyclotron damping playing any role). 
(l.0-exp[-(t/(3.125^ 1 )) 2 ]) factor ensures that 
these driving B± fields ramp up to their maximal values 
in time t — 3. 125a;" 1 . Such driving with B± of 5% of the 
background -Bo excites circularly polarised ion-cyclotron 
(IC) waves, these waves are often misquoted as Alfven 
waves 0, H, 0, [I]. For parallel to B propagation, 
and assuming plasma quasi-neutrality, the relevant 
dispersion relation reads as [ijj : 



k 2 - — 

K R,L — 2 



(UJ =F ^ce)(w ± LU ci ) 



(10) 



where subscripts i?, L refer to the waves with right and 
left circular polarisation, and so are the upper and lower 
signs in the plus and minus. In the frequency range 
lu -C uj c i both left and right polarised waves tend to an 
Alfven wave branch, which satisfies the following disper- 
sion relation [121: 



k 2 -- 



which is the same as 



4:TTn(m e + rrii)c 2 
B 2 



U3 
k 



V A /^l + V 2 /c 2 . 



(11) 



(12) 



Note, that the root in the denominator appears only be- 
cause of the displacement current is kept. At frequencies 
uj ~ 0.3w c i however, the correct term ion-cyclotron wave 
instead of Alfven wave should be used. 

Note that in differ to Refs. 7, 8] we now directly drive 
transverse magnetic fields (B x ,B y ) (which in turn are 
coupled to transverse electric fields (E x ,E y )). Con- 
ventionally, Alfvenic and IC waves are more associated 
with magnetic field oscillation. However, in the kinetic, 
Particle-In-Cell simulation usually variation of electric 
field is used for driving waves, because particles respond 
to electric, rather than magnetic fields (both are natu- 
rally coupled, such as E x and B y represent single Alfven 
wave). In the two- fluid numerical code like the one used 
here, it is possible to use transverse magnetic fields for 
driving of IC wave. 



A. case of mi/m e = 262 

In this subsection we consider case of TOi/m e = 
262.268 rs 262, which is the largest value considered in 
this study. As can be seen in Fig. (3), the generated at 
the left edge (z = 1) IC waves propagate both in the di- 
rections of positive and negative z's. However, because 
of the periodic boundary conditions used (applied on all 
physical quantities) IC wave that travels to the direc- 
tion of negative z's (to the left) re-appears on the right 
edge of the figure. As in all previous phase-mixing sim- 
ulations Alfven velocity is a function of the transverse 
(to the background magnetic field) coordinate, x, i.e. 
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FIG. 3: Contour (intensity) plots of phase-mixed transverse magnetic field B y at times t — 8000/oj ce (left) and t — 16000/aj c 
(right). 
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FIG. 4: Contour (intensity) plots of E z = E\\ at times t = 8000/w ce (left) and t = 16000/w ce (right). 



Va = Va{x) oc l/y/n{x) (see Fig. (2)). Thus as can be 
seen in Fig. (3) the IC wave middle portion travels slower 
than the parts close to the simulation box edge. This cre- 
ates progressively strong transverse gradients and hence 
smaller spatial scales. If resistive effects arc included 
(these are absent here), such configuration usually pro- 
duces greatly enhanced dissipation and IC wave ampli- 
tude decays in space as oc exp(— z 3 ) @, Hj|. The E\\ = E z 
field dynamics is shown in Fig. (4). We gather that E\\ 
is generated only in the regions of density gradients i.e. 
along x = 81 and x = 119 lines. This can be explained 
by analysing right-hand-side (RHS) of Eq.(6). Eu = at 
t = everywhere, however it can only be generated in the 
regions where d x E x ^ 0. The latter is true only in the 
density gradient regions where E x becomes progressively 
oblique propagating. Thus, Eq.(6), derived here for the 
first time, correctly explains the simulated process of En 
generation by IC waves. It should be also noted that 
this equation contains d xx , which correctly accounts for 
the transverse (along x) propagation of the generated En . 
Eu 's longitudinal (along z) propagation due to the mo- 
tion of IC wave along z-axis is indeed corroborated both 



by Fig. (3) and Eq.(7) - note d zz term. Also, note that E\\ 
amplitude at time t = IBOOOlj^ 1 attains value of 3x 1CP 5 . 
This is somewhat smaller value than the one obtained in 
the full kinetic (PIC) simulation This is due to the 

different mass ratios: in the the kinetic (PIC) simulation 
mi/m e — 16, but here it is 262. In dimensional units 
this Eu corresponds to about 0.003 statvolt cm -1 or 90 
V m -1 , i.e. in such electric field electrons would be accel- 
erated to the energy of w 10 keV over the distance of 100 
m. Note, however, that the generated Eu is oscillatory 
in space and time. 

In Fig. (5) we present (V ez — Vi Z ) which is proportional 
to j z , the parallel (electron and ion flow separation) cur- 
rent. It can be seen from this figure that (V ez — Vi 
attains moderate values of w 0.03c. Authors of Ref. 
stated the importance of charge separation before. How- 
ever the cause of En generation is actually electron and 
ion flow separation (see below). The latter is quite dif- 
ferent from the electrostatic effect of charge separation, 
which is inherently a plasma kinetic effect. Electron and 
ion flow separation is a fluid- like (non- kinetic) effect, be- 
cause their distribution functions remain Maxwellian at 
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FIG. 5: Contour (intensity) plots of (Vez — Viz) oc j z at times t = 8000/cJ ce (left) and t — 16000/w ce (right). 



all times. 



B. parametric study for different mj/m e 's 

In this subsection we perform parametric study with 
different mass ratios mi/m e . This is particularly use- 
ful for understanding the physics of parallel electric field 
generation. This is needed because performing realistic 
mass ratio numerical simulation, mi/m e — 1836, is com- 
putationally challenging. 

Numerically challenging issues arise from the increase 
in mass ratio are the following: 

(i) because our intention was to use a driving IC wave 
with uifi = 0.3uj C i, when changing mass ratio, Wd oc 1/toj. 
Thus increasing mass of the ion leads to decrease in driv- 
ing frequency. Since in all our numerical simulations we 
intended to have final simulation time of 3 driving IC 
wave periods, i.e. tfi na i = 3 x (2ir/ujd) oc m t . Thus, 
in turn, increase in ion mass leads to increase in final 
simulation time. E.g. for the cases rrii/m e = 45.9, 91.8, 
183.6, 262.286 considered here, the final simulation times 
were 2884, 5768, 11536 and 16000 respectively, i.e. 

i final OC 77lj. 

(ii) Increase in simulation time leads to the compulsory 
increase in the simulation domain size, L, in the direction 
along the magnetic field (direction of IC wave travel), i.e. 
L = 2 x (Va/c) x {it f^ai) oc (1/^/toT) x m l oc ^/ml 

Thus, as far as the slowdown of the numerical code 
is concerned, the combined slowdown effect of the above 
two factors is rni^/mi. 

In Fig. (6) we present time evolution of the amplitude 
of the generated parallel electric field, which we define 
as max(\E z (x, z, t)\) 7 i.e. at every time step we choose 
one point over whole simulation domain at which mod- 
ulus of parallel electric field is the largest. This is, as 
it was shown above, along the strongest gradient lines 
x = 81J and x — 119(5. It can be seen from the graph 
that: (i) level attained by parallel electric field amplitude 
decreases with the increase in the mass ratio; (ii) rate at 
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t 

FIG. 6: Time evolution of the amplitude of the generated 
parallel electric field, defined as max(|_B z (a;, z, t)\). The thick- 
est (upper left) line corresponds to m;/m e = 45.9; while lines 
with decreasing thickness correspond to 91.8, 183.6, 262.286 
respectively. 



which the final amplitude is reached (the averaged slope, 
essentially) also decreases with the increase in the mass 
ratio. 

Exact behaviour of the final attained parallel electric 
field amplitude (within 3 periods of the driving ion cy- 
clotron wave) as a function of mass ratio is shown in 
Fig. (7). Essentially this is a plot of the last data points 
(which are four) in the previous Fig. (6) as a function 
of m,/m e , i.e. E* = max(\E z (x, z, tf ina i)\) vs. mi/m e . 
The dashed line corresponds to the fit 0.0085/(mi/m e ). 
Plotting such graph is very useful to establish the trend. 
Interestingly we see that that amplitude attained by E\\ 
decreases linearly with inverse of the mass ratio rrii/m e . 
The x-range in Fig. (7) is mi/m e = 30 — 1836, so that 
rightmost point of the dashed line enables us to grasp En 
for the case of realistic mass ratio (i.e. 1836). We thus 
gather that En = 0.0085/1836 = 4.630 x 10~ 6 which is 
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FIG. 7: The final attained parallel electric field amplitude 
(generated within 3 periods of the driving ion cyclotron wave) 
as a function of mass ratio. Data points correspond to the 
last 4 data points in the previous Fig. (6). The dashed line 
corresponds to the fit 0.0085/(m;/m,;). This is a log-log plot. 



0.035 




FIG. 8: Time evolution of the amplitudes of the generated 
parallel flows of electron fluid, defined as max(| V e z(x, z, t)\). 
The thickest (upper left) line corresponds to m,i/m e = 45.9; 
while lines with decreasing thickness correspond to 91.8, 
183.6, 262.286 respectively. 



4.630 xl(T 6 x E a = 4.7 x 1(T 4 statvolt cm" 1 or 14 Vm" 1 . 

In Figs. (8) and (9) we plot the amplitudes of the 
generated parallel flows of electron and ion fluids which 
we define as max(|V ez (a;, z, t)\) and max.{\Vi Z (x, z, t)\) re- 
spectively, i.e. at every time step we choose one point 
over whole simulation domain at which moduli of paral- 
lel flows of electrons and ions are the largest. Again, this 
occurs somewhere along the strongest gradient lines x = 
815 and x = 1195, because parallel electron and ion flows 
are confined to the strongest density gradient regions. 
Four different lies in each figure show the cases rrii / m e — 
45.9, 91.8, 183.6, 262.286 considered. We gather from 
Figs. (8) and (9) that an increase in mass ratio does not 
have any effect on final parallel (magnetic field aligned) 
speed attained by electrons. However, parallel ion ve- 
locity decreases linearly with inverse of the mass ra- 
tio rrii/m e . To investigate this more quantitatively, in 
Fig. (10) we plot the ratio of final attained electron and 
ion flow amplitudes (within 3 periods of the driving ion 
cyclotron wave) as a function of mass ratio. Essentially 
this is a plot of the ratio of last data points (which are 
four) in the Figs. (8) and (9) as a function of mi/m e , i.e. 
Vratio = xaax.(\V ez (x,z,tfi na i)\)/max(\V iz (x,z,tfi n ai)\) 
vs. mi/m e . The dashed line corresponds to the fit which 
is (mi/m e ) with a slope of 1. i.e. parallel velocity ratio 
of electrons and ions scales directly as mi/m e . 



III. DISCUSSION 

The conclusions that follow from the collective analysis 
of Figs.(6)-(10) initially may seem counterintuitive. On 
one hand maximal attained Eu amplitudes drop off as 
l/rrii (Figs.(6)-(7)). On the other hand, electron flow 
maximal attained amplitudes do not depend on rrii (they 
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FIG. 9: The same as in Fig. (8), but for ions. 



all are circa 0.03c, see Fig. (8)), while ion flow maximal 
attained amplitudes drop off as l/rrii (Figs. (9)-(10)). 
Thus one might expect that more massive ions should 
produce a bigger E\\ (since separation of electron and ion 
fluids is the source of Eu and that separation is expected 
to be largest in the case of more massive ions, as they 
are slower). In fact, this is what would be expected if 
the polarisation drift produced by the driving IC wave is 
the cause of parallel electric field generation [5, 6]. The 
latter two references use the following polarisation drift 
current: 



3x 



miUi dE± 
~B 2 dt 



(13) 



where symbols have their usual meaning. The latter 
equation implies that Eu then should increase with oc 
rrii. However, in Figs.(6)-(7) we see completely oppo- 
site E\\ tx \ I mi scaling. These results can be interpreted 
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FIG. 10: The ratio of the final attained electron and ion 
flow amplitudes (within 3 periods of the driving ion cy- 
clotron wave) as a function of mass ratio, i.e. Vratio = 
ma,-x(\Vez(x,z,tf ina i)\)/ma,x(\Vi z (x,z,tfi„ a i)\) vs. mi/m e . 
The dashed line corresponds to the fit which has a slope of 1 
(note the direct correlation). This is a log- log plot. 

(reconciled) as following: (i) ion dynamics plays no role 
in the En generation, i.e. polarisation drift has no effect 
in contrary to the claims of Refs. [H, Q; (ii) decrease in 
the generated parallel electric field amplitude with the in- 
crease of the mass ratio /m e is caused by the fact that 
uod = 0.3w C i oc l/rrn is decreasing, and hence the electron 
fluid can effectively "short-circuit" (recombine with) the 
slowly oscillating ions, hence producing smaller En which 
also scales exactly as l/m*. 

In summary, indeed, electron and ion fluid separation 
is causing En generation, but polarisation drift current 
produced by the driving IC wave plays no role. 

Interestingly, by comparing Figs. (8) and (9) we gather 
that electron fluid is efficiently accelerated by the gener- 
ated En to the velocities of up to 0.03c, while ion fluid 
due to its larger inertia is much less mobile (0.0001c — 
0.0007c). This confirms yet another conclusion that was 
made in Refs. [7], [8( which employed full kinetic simula- 
tion. 

It should be noted that since here we use two-fluid 
approach the generated En cannot change the distribu- 
tion function, which obviously remains Maxwellian, while 
in the previous kinetic simulation of a similar system it 
produced bumps in the distribution function as the elec- 
trons residing on the magnetic field lines with the density 
gradients get efficiently accelerated (see e.g. Fig. (4) in 
Ref.l). 

IV. SUMMARY 

We studied the generation of parallel electric fields by 
means of propagation of IC waves in the plasma with the 
transverse density inhomogeneity. By using simpler, than 
kinetic [EH SB], two-fluid, cold plasma linearised equa- 



tions, we show for the first time that En generation can 
be understood by an analytic equation that couples En to 
the transverse electric field. It should be noted that the 
generation of En is a generic feature of plasmas with the 
transverse density inhomogeneity and in a different con- 
text this was known for decades in the laboratory plasmas 
(15 , [l4| . We prove that the minimal model required to 
reproduce the previous kinetic results of E\\ generation 
is the two-fluid, cold plasma approximation in the linear 
regime. In the latter, the generated En amplitude attains 
values of 14 Vm _1 for plausible solar coronal parameters 
and realistic mass ratio of mi/m e = 1836. By consider- 
ing the numerical solutions for (V ez — Vi Z ), we have shown 
that the cause of E\\ appearance is the electron and ion 
flow separation (which is not the same as electrostatic 
charge separation) induced by the transverse density in- 
homogeneity. 

We also investigate how E\\ generation is affected by 
the mass ratio and found that amplitude attained by En 
(within 3 periods of the driving ion cyclotron wave) de- 
creases linearly as oc 1/toj. This result contradicts to 
the earlier suggestion by Genot et al (1999, 2004) that 
the cause of E\\ generation is the polarisation drift of 
the driving wave, which would suggest E\\ oc rrij scaling. 
Increase in mass ratio does not affect the final parallel 
(magnetic field aligned) speed attained by electron fluid. 
However, parallel ion velocity decreases linearly as 1/rrii, 
this means that the parallel velocity ratio of electrons 
and ions scales directly as m*. 

It should be noted that when plasma density is homo- 
geneous, no E\\ generation takes place; and this is cor- 
roborated both by numerical simulations (not presented 
here) and agrees with the Eq.(6) (when n — const, the 
RHS of Eq.(6) is zero at all times as E x and B y do not 
propagate obliquely). Our model also correctly repro- 
duces the previous kinetic results [E IE 0> HI that only 
electrons are accelerated (along the background magnetic 
field), while ions practically show no acceleration. 



V. APPENDIX 

Animations 1,2, and 3 show movies corresponding to 
Figs. (3), (4) and (5) respectively. Note that horizontal 
axis indicates 500 grids, while the real simulation value 
is 2500. This is simply to reduce movie size, i.e. every 
5-th point along the field was included in these movies. 
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Note added in proof: The typical values of the Dre- 
icer electric field on the corona is a few x 10 3 V m" 1 [ill ], 
which implies the obtained E\\ in our model exceeds the 
Dreicer value by at least four orders of magnitude, per- 
haps enabling the electron run away regime. This would 
imply that our model is more relevant to the accelera- 
tion of solar wind, rather than solving the coronal heat- 



ing problem. Essentially acceleration of electrons would 
dominate over the heating as such. However, this seems 
uncertain because electron and ion fluid separation can- 
not go on (build up) forever, and some sort of discharge 
should eventually take place. At any rate, similar kinetic 
simulations have shown [|| (see their figure 11) thatwave 
energy is converted into particle energy on timescales of 
lO 3 ^^ 1 (mind that the latter number is likely to be de- 
pendent on the mass ratio mi/m e ). 
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